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L'invention conceme notamment un precede d'identification 
autodidacte d'un nombre de sources P potentiellement superleur ou egal au 
nombre N de capteurs de Tantenne de reception. 

EUe peut etre utilisee par exemple dans un contexte multi- 
emissions a Bande etroite. 

Elle est utilisee par exemple dans un reseau de communications, 

EUe trouve son application notariiment dans le domaine des 
radiocommunications, des telecommunications spatiales ou de Tecoute 
passive de ces liaisons, dans des gammes allant de la VLF a la EHF. 

La figure 1 schematise un exemple de reseau comprenant 
plusleurs capteurs de reception, chaque capteur recevant un ou plusieurs 
emetteurs de radiocommunication, de directions d'arrivee differentes. t 

Chaque capteur regoit ime source avec une phase et ime 
amplitude dependant de Tangle d'incidence de celle-ci et de la position du 
capteur. La figure 2 schematise un exemple de parametrage de la direction 
dune source. Cette direction est parametree par 2 angles correspondant 
aux angles d'azimut 6 et d'elevation A, 

Depuis maintenant pres d*une quinzaine d'annees, de 
nombreuses techniques d'identification autodidacte de signatures ou 
vecteurs directeurs de sources, supposees statistiquement independantes, 
ont ete developpees sous Thypothese d'un nombre de sources P inferieur ou 
egal au nombre de capteurs N, decrites dans les references (1][3]I7I citees en 
fin de la description. Toutefois pour de nombreuses applications pratiques 
telles que les radiocommunications HF, le nombre de sources regues par les 
capteurs augmente en particulier avec la largeur de bande des recepteurs et 
le nombre de sources P .peut des lors devenir superieur au nombre de 
capteurs N. Les melanges associes des sources sont alors qualifies de sous- 
determines. 

Un certain nombre de methodes dldentification autodidacte de 



1 er depot 




2 

melanges sous-determines de sources^ a bande etroite pour le reseau ont ete 
developpees tout recemment et decrites dans les references [2] [7-8] et [10]. 
Les methodes proposees dans les references [2] et [7-8] e>q)loltent 
Tinformation contenue dans les statistiques d'ordte-.4 (FO : Fourth Ordpr) 
5 des slgnaux regus sur les capteurs tandis que la methode proposee dans la 
reference [10] exploite Tinformation contenue dans une des fonctions 
caracteristiques des signaux regus. Toutefois, ces methodes presentent de 
grosses limitations dans la perspective d'une mise en oeuvre operationnelle. 
En effet, la methode de la reference [2] est tres difiQcile a mettre en oeuvre et 

10 n'assure pas Fidentification des sources ayant le meme kurtosis (cumulant 
normalise d'ordre 4). Les methodes des references [7-8] supposent que les 
sources sont non circulaires et conduisent a des resultats peu fiables en 
pratique. Finalement la methode de la reference [10] a ete developpee 
xmiquement pour des melanges de sources a valeurs reelles (non complexe). 

15 L'objet de la presente invention conceme notamment un nouveau 

procede dldentification autodidacte d*un melange sous-determine de 
sources a bande etroite pour le reseau. Le procede permet notamment 
dldentiiier jusqu'a - N + 1 sources a partir de N capteurs identiques et 
jusqu*a N2 sources avec N capteurs differents, en supposant uniquement 

20 que les sources ont des trispectres differents et des kurtosis non nuls et de 
meme signe (cette demiere hypothese est pratiquement toujours veriflee 
dans le contexte des radiocommunications). 

L*invention conceme un procede d'identification autodidacte a 
25 Tordre 4 de signatures d*au moins deux sources dans un systeme 
comportant un nombre de sources P et un nombre N de capteurs de 
reception recevant les observations, lesdites sources ayant des trl-spectres 
differents. U est caracterise en ce qu'il comporte au moins les etapes 
suivantes : 

30 • une etape de blanchiment a Tordre 4 des observations regues sur les 
capteurs de reception afln d*orthonormaliser les vecteurs directeurs des 
sources dans les matrices de quadricovariance des observations 
exploitees. 
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• line etape de diagonalisation'^ conjointe de plusieurs matrices de 
quadiicovariance blanchles afin d'identifier les signatures spatiales des 
sources. 

Le nomb^e-'de sources P est par exemple superieur au nombre de 
capteurs N. 

Le procede peut etre utilise dans un reseau de communication. 
Le procede selon Tlnvention presente notamment les avantages 

suivants : 

• n permet d'identifler un nombre de sources superieur au nombre de 
capteurs : 

• pour des capteurs identiques : - N + 1 sources ayant des 
trlspectres differents. des kurtosis non nuls et de meme signe» 

• pour des capteurs differents {reseau a diversite de polarisation et/ou 
a diversite de diagrarmne et/ou avec du couplage, etc..) sources- 
ayant des trispectres differents, des kurtosis non nuls et de meme 
signe, 

• n est robuste a du bruit Gaussien meme correle spatialement» 

• II permet une goniometrie de chaque source identifiee, a partir d'un 
modele de front d'onde plaque sur la signature, avec une resolution 
potentiellement superieure aux methodes existantes, 

• II permet ridentiflcation de I (N^-N+l) sources cyclostationnaires si les 
capteurs sont identiques et Ix N2 sources cyclostationnaires si les 
capteurs sont differents : a diversite de polarisation et/ou a diversite de 
diagramme et/ou avec du couplage. ou I est le nombre de frequences 
cycliques traitees, 

Au moyen d'un critere de performances, il permet d6valuer 
quantitativement la qualite d'estimation du vecteur directeur de chaque 
source, tme comparaison quantitative de deux methodes pour 
Tidentification d'une source donnee, 

A Taide d'une etape de selection des frequences cycliques, il permet de 
tralter un nombre de sources superieur au nombre de sources traitees 
par la methode de base. 
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D'autres caracterlstiques et avantages de rinvention apparaitront 
mieux a la lecture de la description annexee des figures qui representent : 

• La figure 1 un exemple de reseau de communications, 

• La figure 2 une representation des parametres'd^iirie source, 

• La figure 3 un schema fonctionnel du procede selon Tinvention, 

• La figure 4 un exemple de flltrage spatial, 

• Les figures 5 et 6 des exemples de variation du critere de performance 
en fonction du nombre d'echantillons observes comparant les 
performances du procede a 2 methodes de Tart anterieur. 

• Les figures 7 a 9 trois variantes de realisation du procede decrit a la 
figure 3 mettant en oeuvre une selection des frequences cycliques. 

Afin de mieux faire comprendre Tobjet de rinvention, Texemple 
qui suit est donne a titre illustratif et nullement limitatif pour un reseau de 
radiocommunication en contexte de multiemissions a Bande etroite, avec 
des sources ayant des tri-spectres (de cumulants) differents. 

Chaque capteur du reseau, compose de N capteurs, regoit un 
melange de P sources a Bande Etroite (BE) supposees statistiquement 
independantes. Sous cette hypothese, le vecteur des enveloppes complexes 
des signaux en sortie des capteurs s'ecrit : 

p 

x(t) =2 SpW Up + b[t) = A sit) + h(i) (1) 

ou Sp(i9 est le signal de la p^me source ainsi que la peme composante du 
vecteur s(t ), h(t ) est le vecteur bruit suppose gaussien de covariance 
quelconque, Op est la signature ou le vecteur directeur de la p^^^ source et A 
(N X P) est la matrice des vecteurs ap (vecteurs directeurs dies sources). 

Llnvention a notainment pour objet d'identifier les vecteurs 
directeurs Up de chacune des sources lorsqu'en particulier le nombre de 
sources P est potentiellement superieur au nombre de capteurs iV. 

A partir de cette identification, il est ensuite possible d'appliquer 
des techniques d^extraction des sources par filtrage spatial des 
o.bservations. ^extraction autodidacte a notamment pour objectif de 
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restituer les . signaux informationnels vehicules par les sources en 
n'exploitant (en fonctionnement normal) aucune Information a priori sur 
celles-ci. 

Statistique d'ordre^4 

5 Le procede selon Tinvention exploite les statistiques d'ordre 4 des 

observations correspondant aux moyennes temporelles. 0x(ti,T2»T3) 
=<0;c(ti,T2,T3)(O>, sur un horizon d'observation infini, de certalnes matrices 
de quadrlcovariance. 0,v(ti.T2.T3)(O. de dimension [N^xN^]. Les elements, 
Qx(xi,T2.T3)It J. fc. ZKU. de ces matrices sont par exemple definis par la 
10 relation : 

Qx(Ti,X2,X3)[i,j\ K (][(] = CumCY/(0, xj{i -x,)*, x^Ht -Tj)*, -T3)) ' (2) 

ou * est le symbole complexe conjugue. Xitfl est la i composante du 
vecteur , <.> est I'operation de moyennage temporel sur un horizon 
15 d'observation infini et (ti,T2.T3) est un triplet de retards. En supposant que 
Qx{Tui2.X3)[l J, /c, g est Telement [N(i -1) + j, N{k - 1) + Z] de la matrice 
^ 0,v(ti,T2,T3), en supposant que le bruit est gaussien et en utilisant 
Texpression (1) dans Texpression (2), la matrice de quadrlcovariance 
9x(ti,T2»T3) s*ecrit de la fagon suivante : 

20 - Q.{Ti, T2, Tj) = (A ®A*) Qs(Tj. T2. Ts) {A ®Af (3) 

ou 0s{ti,T2.T3) est la matrice de quadrlcovariance rtloyennee de s(<J de 
dimension (P^xF^], A=l ai ... ap], ® est le produit de PCronecker et h designe 
le transpose et conjugue. Sous Thypothese de sources statistiquement 
independantes la matrice Q5(ti,T2,T3) est composee d'au moins - P zeros 

25 et Texpression (3) se simplifie de la maniere suivante : 

p 

Qx(Xl ,T2,T3) = 2 ''^2,X3) (a^® Op *) {Op ® Up (4a) 

= Aq C,(Ti,T2,T3) Aq " (4b) 
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ou Aq est line matrice de dimension {N^x Pi definie , par Aq = [(ai®ai*), 
{ap®ap*)], Cs{ti,T2.T3) est une matrice diagonale de dimension {P x F5 definie 
par Cs(ti,T2,T3) = diag[ci(xi,X2,T3)...., CpCxuTa.xs)] et 6u^iS,(ti,X2.X3) est defini par : 

Cp(Xi,X2,X3) = <Cum(5p(0, Sp{t'X{)\ Sp{t'X2)% Sp(t-Xi))> (5) 

5 L'expression (4b) a ime structure algebrique similaire a celle de la 

matrice de correlation des observations utilisee dans Talgorithme SOBI 
decrit dans la reference [1]. On notera dans la suite Qx = 0x(O, 0, 0), Cp = 
Cp(0, 0. 0), Cs = Cs(0, 0, 0) pour en deduire de la relation (4b) : 

0x= Aq Cs ApH (6) 

^0 On suppose, dans la suite, que le nombre de sources P est tel 

que P < JV^, que la matrice Aq est de rang plein, que les cumulants 
moyennes Cp , 1 < p < P, sont non nuls (sources non gaussienne) et de 
meme signe et que pour tout couple iUJ) de sources il existe au moins un 
triplet de retards (xi,X2,X3) tel que | Xi | + 1 X2 1 + 1 xs | 5^ 0 et 

15 Ci(Xi,X2.X3) / \Ci\ ^ Cj(Xi,X2,X3) / \Cj\ (7) 

Etape de blanchiment a I'ordre 4 

La premiere etape du procede selon Tinvention, denomme 
20 FOBIUM, consiste a orthonormaliser, dans la matrice de quadricovariance 
Qx de Texpression (6). les colonnes de la matrice Aq, considerees comme les 
vecteurs directeurs vlrtuels des sources pour le reseau de capteurs 
considere. Dans ce but, le procede considere la decomposition en elements 
propres de la matrice hermitienne Qx, de rang P, donnee par 
25 a = . (8) 

ou Ax est la matrice diagonale reelle, de dimension (P x JFO, des P valeurs 
propres non nuUes de Qx. et Ex est la matrice de dimension [N^xF) des 
vecteurs propres prthonormes associes. Pour une matrice Aq de rang plein, 
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on peut montrer qu'il y a equivalence entre supposer que les kurtosls des 

sources sont de meme signe e (e = ±1) et supposer que les valeurs propres 

de Ax sont aussl.<^e..menie signe e. Dans ce contexte. il est possible de 

construire la matrice de blanchiment suivante Tde dimension (P x iV-^): 
5 T-iAS'^'E.'' (9) 

La matrice de blanchiment de dimension (P x AT^) est definie a 
partir de la racine cairee de la matrice diagonale reelle de dimension (P x P) 
des P valeurs propres non niUles de la matrice de quadricovariance et du 
transpose de la matrice des vecteurs propres associes de dimension (P x N^) 

10 ou (Aa) est Tinverse de la racine carre de A^. Des expressions (6) et (8) on 

en deduit que : 

sTQ^f = TAoi£Cs)AQ''f' - Ip (10) 
ou Ip est la matrice identite de dimension (P x et ou eCs = diag[|ci|,.... |Cp|]. 
Cette demiere expression montre que la matrice TAq {eCsV^^ de dimension (Px 
15 P) est une matrice unitaire U. On en deduit alors que : 

-1/2 

TAq = UieCs) (11) 

Etape d'identification a Fordre 4 

On deduit des expressions (4b) et (1 1) que : 

20 WiXuX2M=-Ta(ruX2,X3)T^==U(£Cs)'^^^^ {eC,)'^^ if (12) 

ou XV{ti.T2.T3) est la matrice de quadricovariance blanchie a Tordre 4 par la 
matrice Q^. Cette deitiiere expression montre que la matrice unitaire U 
diagonalise les matrices TQAti,T2.x^) T» et que les valeurs propres associees 
sont les termes diagonauLX de la matrice diagonale (eCs)^^^ Cs(ti,X2,T3) (eCs)'^^. 
25 Pour un triplet de retards donne (Ti,t2.T3), la matrice U est linique a une 
permutation pres et a une matrice diagonale unitaire pres, lorsque les 

—1/'' — i/'' 

elements de la matrice (eCs) CsixuX^M (eCs)"^sont tous differents. Dans le 
cas contraire, le procede utilise un ensemble de K triplets (ti^T2^T3'^) , 1< fc < 
K, defini de la maniere suivante : pour tous couples de sources (i, jO. il existe 
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au moins un triplet (xi^xa'^.Ta'^). tel que la conditioh de Tequation (7) soit 
verifiee. Dans ces conditions, la matrice unitaire U est la seule matrlce Usoi 
qui, a une permutation et une matrice diagonale-linitalre pres, diagonalise 
conjointement les K matrices TkQJtiKt2Kt3^)xT^. En consequence, la matrice 
5 Usoi. solution du probleme precedent, s'ecrit en fonction de la matrice 
unitaire U de la maniere suivante : 

Usoi =U ATI (13) 

ou A et n sont respectivement la matrice diagonale unitaire et la matrice de 
permutation citees precedemment. A partir des equations (11) et (13), il est 
0 possible de deduire la matrice Aq a une matrice diagonale unitaire et de 
permutation pres. qui s'exprime par : 

7^ f/soi aJ"^ f/soi (sCs)^^ A n (14) 

ou 'f est la pseudo-inverse de la matrice T. Chaque colonne > (1< f < f), de 
la matrice 'f Usoi correspond a Tun des vecteurs jjiq |Cq|l/2 (aq®aq*), 1< q < P, 
ou liq est un scalaire complexe tel que lml = 1. En consequence, en 
transformant chaque colonne bi de la matrice 7^ Usoi en une matrice Bi de 
dimension (iVx W) tel que Bilij] = b?((i-l)JV+ j) (1 < u j < iV) on en deduit que: 

1/2 J, 

Bi^]i^\Cg\ a,^ag pour (l<l,q<P) (15) 
La matrice Bi est construite a partir du vecteur bi et depend d'un 
scalaire complexe, de la raclne carree du cumulant et du vecteur directeur 
de la q^eme source et de son conjugue. 

Dans ce contexte le vecteur directeur a<j de la qf^^^^ source est associe au 
vecteur propre de Bi associe a la plus forte valeur propre. 
Resume du principe de I'lnvention 

En resume, les differentes e tapes du precede selon I'invention 
comportent au moins les etapes suivantes : 

pour L observations vectorielles regue au cours du temps : x[lTe) {l< l< L). 
ou Te est la periode d'echantillonnage. 

Estimation 
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Etape 1: Estimer, par la matrice de quadricovariance Qx. a partir des L 
observations x{lTe), en utilisant un estimateur non bialse et 
asymptotiquementjconsistant. Suivant la nature des sources. Testimateur 
est adapte comme suit : 
5 • Cos stationnaire et centre : Estimateur empirique utilise dans la reference 
[3]. 

• Cos cyclo-stationnaire et centre : Estimateur mis en oeuvre dans la 
reference [10]. 

• Cos cyclO'Stationnaire et non-centre : Estimateur mis en oeuvre dans la 
10 reference [11]. 

Blanchiment 

Etape 2: Decomposer en elements propres la rnatrice de quadricovariance 
estimee estimer le nombre de sources P et restreindre cette 

decomposition en elements propres. aux P composantes principales : tx 
15 Ax kjf^r ou kx est la matrice diagonale contenant les P valeurs propres de 
plus fort module et kx est la matrice contenant les vecteurs propres 
associes. 

Etape 3 : Construire la matrice de blanchiment : f= (Aj'^'^ jfe^. 
20 Selection des triplets 

Etape 4 : Selectionner K triplets de retards [ii\x2^.i^^) ou 

|Tl'^| + |T2^| + |T3'^|9tO. 

Estimation 

Etape 5 : Estimer. par ^.v(ti^T2^T3'^), K matrices de quadricovariance Qx 
25 (ti^T2^T3'0. Conrnie dans Tetape 1 cette estimation depend notamment des 
hypotheses faites sur les observations: . 

• Cos stationnaire et centre : Estimateur empirique utilise dans la reference 
13], 



1 er depot • 



10 

• Cos cyclo-stationnaire et centre : Estimateur mis en oeuvre dans la 
reference [10]. 

• Cos cyclo-stationnaire et non-centre : Estirtiateur^rnis en oeuvre dans la 
reference [11]. 

5 Identification 

Etape 6 : Calculer les matrices t ^Jzi^.'i2Kt3^) 1^ et estimer. par fJsoi, la 
matrice unitaire Usoi par diagonalisation conjointe des K matrices 

Etape 7 : Calculer f 'tlsoi=[hi...hp] et construire les matrices hi de 
0 dimension (iVx N). 

Etape 8 : Estimer, par 2Lp, les signatures a^ {l< q < Pi des P sources en 
appliquant ime decomposition en element sur chaque matrice hi„ 
Applications 

A Tissue de Fetape 8, le precede a identifle les vecteurs directeurs 
5 de P sources non gaussiennes ayant des tri-spectres diiferents avec des 
kurtosis de memes signes. P<N2 et P peuvent atteindre - N + 1 ou 
selon le t5rpe de capteurs utilises. 

A Taide de ces informations, le procede peut mettre en oeuvre une 
methode de goniometrie ou im filtrage spatial d'antennes. 

Une methode de goniometrie permet de determiner la direction 
d'arrivee des sources et plus precisement les angles d'azimut, Gm pour la 
goniometrie en ID et azimut et site (Gm, Am) pour la goniometrie 2D. 

La figure 4 represente un filtrage spatial d'antennes, pour des 
structures spatiales de filtrage. II permet notamment d'optimiser la 
reception d'une ou de toutes les sources presentes par filtrage spatial des 
observations. Lorsque plusieurs sources sont d'interet pour le recepteur, on 
parle de techniques de separation de sources. Lorsque qu'aucune 
information a priori n'est exploitee sur les sources, on parle de techniques 



1 er depot 



11 



autodidactes. 

Verification de la qualite des estimees 

Selon tine variante de realisation, le precede comporte line etape 
5 permettant revaluation quantitative, pour chaque source, de la qualite de 
Tidentification du vecteur directeur associe. 

Ce nouveau critere permet la comparaison intrinseque de deux 
methodes d'identificaUon pour la restitution de la signature d'une source 
particuliere. Ce critere, pour le probleme de ndentificatlon, est Textension 
10 de celui propose dans (5] pour r extraction. II est defini par le P-uplet : 

D(A, A) = (ai, a2, , a^) (16) 

ou 

ot^ = min [d(ap, fl/)] - (17) 
et ou d(u,ii) est la pseudo-distance entre les vecteurs u et u, tel que : 



15 



u V 



Dans les simulations des figures 5 et 6. on est en presence de 
P=6 soxirces statistiquement independantes regues sur im reseau circulaire 
de 7V=3 capteurs de rayon r tel que r/A=0.55 (A: longueur d'onde). Us 6 
sources sont des sources QPSK non filtrees ayant 20dB de rapport signal 
20 sur bruit avec une periode symbole T = 4Te. ou Te est la perlode 
d'echantillormage . 

Les incidences des sources sont telles que &i=2.16°, 4=25,2°, 
6b=50°. a=272.16°. ^=315.36°. 6fe=336.96° et les frequences porteuses 
associ^es verifient A/i Te=0. A/2T«=l/2. A/3 Te-1/3, A/4 Te=l/5, A/5 Te=l/7 et 
25 A/e Tc=l/ll. Les methodes JADE (3), SOBI [1] et FOBIUM selon Tinvention 
sont appliques et les performances ou, pour q=1...6 sont evaluees apres un 
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moyennage sur 1000 realisations. Pour la methode FOBIUM on choisit K=4 
triplets de retards (ri^, 7-2^^3'^) ou ri'^^/cTeet T2^-r3^=^0. 

Sous les hypotheses precedentes. la figure 5 montre la variation 
5 de a2 (performances de la 2ieme source) en sortie des separateurs JADE (b), 
SOBI (c) et FOBIUM (a) en fonction du nombre L d'echantUlons. Les courbes 
montrent d'une part que les methodes JADE et SOBI ont des difflcultes a 
identifier le vecteur directeur de la 2^^™* source dans un contexte de 
melange sous-determine et que d*autre part la methode FOBIUM a de trts 
10 bonnes performances. 

La Figure 6 montre, dans le meme contexte, les variations de 
tous les ap (1 < p < 6) en sortie de la methode FOBIUM en fonction de L. La 
coin-be (indice p) est associe a la p'^e source. On remarque que tous les 
15 coefficients ap convergent vers zeros et que asymptotiquement les vecteurs 
directeurs sont parfaitement identifies. 

Variantes de realisation cyclique 

Les figures 7 et 8 representent deux exemples de variantes de 
20 realisation selon I'invention denommees FOBIUM cyclique. 

L'idee consiste notamment a introduire une selectivity par les 
ft-equences cycliques dans le procede presente ci-dessus et vise en 
particulier a identifier de maniere autodidacte et avec une capacite de 
traltement plus importante, des melanges sous-determines de sources 
25 cyclostationnaires. 

La difference majeiire entre les etapes 1 a 8 expllcitees ci-dessus 
et cette vaxlante de realisation est la mise en ceuvre d'lme etape dlsolation 
cyclique des sources par discrimination selon leurs frequences cycliques a 
Tordre 4. De cette fagon, H est possible d'identifier separement les sources 
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associees a un meme paramdtre 'cyclique k Tordre 4 sans etre perturbe par 
les autres sources traitees separement. 

Les deixx variantes presentees aux figures 7 et 8 peuvent etre 
mises en oeuvre en reiterant le processus d'isolation cyclique sur les 
5 « autres sources avec d'autres parametres cycliques. Le processus 
d'isolation cyclique peut s'appliquer plusieurs fois dans une troisieme 
version illustre sur la figure 9. 

Statistiques c^cliques d'ordre 4 

10 Les statistiques d'ordre 4 cycliques des observations ou des 

signaux capteurs utilisees sont caracterisees par les matrices de 
quadricovariance cyclique 0x^(a,Ti,T2,T3) de dimension (JV^xiV^) ou les 

elements 0x^(a.TuT2.T3)[i, j. fc, i] , sont deflnis par : 
aka,Ti,T2,T3) [ij, /] = <Cum(x/(0. xj(t 'Xx)\ xkit -T2)*, xi(t -T3)) exp(-j27ta/)> 
15 fir^(a,Ti,X2,T3) [ij, I\ ^ <Cum(x/(0, xj{t -T,) * , Xkit -T2) , xi{t -T3)) exp(-j27ca/)> 

a\a,T,,T2,T3) [ij, KI] ^ <Cum(xX/), xj(j -T,) , xi^t -T2) , xi(t -T3)) exp(-j27:a/)> ( 1 9) 

On note que 2rlCa,Ti,X2,X3) est associe au f "^""^ moment d'ordre 4. En ppsant 
que fir^(a,T,,T2,T3)[/, y, I] est Telement [A^(/ -1) + 7, A^(^ - 1) + q de la matrice 
gr^(a,T,,T2,T3) ct supposant que le bruit est gaussien la matrice 2r^(a,T|,T2,T3), 
20 s'ecrit de la fagon suivante en utiUsant (1) et (19) : 

Qx\cL,Ti, Ts) = (A ®^') Qs\a,Tj, T2. Ts) (A 0AY 

2x^(cx, r;, r^, Ts) = ) 0/(a, zy, ^2, rj) (A ®Af 

Qx\a,Tj, T2. Ts) = (A®A) e/(cc, z), ^2, Ts) (A ©Af (20) 

Ou 0s^(a.Ti,T2,T3) est une matrice de quadricovariance cyclique de s{Q de 
25 dimension (P^xi^), ®est le produit de Kronecker et designe le transpose. 
Sous rh3rpothese de sources statlstlquement independantes la matrice 
0s^a,Ti,T2.T3) est composee d'au moins - P zeros et Texpression (20) se 
simplifie de la maniere suivante : 
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r ■ ^ ■ . 

P 

eAot,Ti,T2,T3) = 2] c/(a,T,,T2,X3)(«/,®a;,)(a;,®a;,f = 5^ C,\a,T,,T2,T3) 5/ (21) 

ou et Bq sont des matrices de dimension (N^ x P ) definie par Aq = 

5 [(ai®ai*) [ap®ap)] et Bq = [(ai®ai) , (ap®ap)], Cs'"(a,Ti,T2,T3) est une 

matrice diagonale de dimension (P x definie par Cs^(a,Ti,T2,T3) = 

diagIci^(a,Ti,T2,T3) Cp^(a,xi,X2,T3)l et ou Cp^(a,xi,X2,X3) est defini par : 

Cp\a ,x,,X2,X3) = <Cum{sp{t\ sp{t-Xx)\ sp(t-X2)\ Spit-x^)) exp(-j27icx/)> 
Cp\a ,x,,X2,X3) = <Cum(spO\ Sp(f-Xx)\ Spit-Xj) , Sp{t-x^)) exp(-j27ccx/)> 
1 0 Cp\cL ,Xi,T2,X3) = <Cum(5p(/), sp{t-x{), sp(t-X2) , spit'X^)) exp(-j2Tca/)> (22) 

On note que la quadricovariance classique de (6) verifie aussi que 
Qx = Qx (0, 0, 0, 0). Cp = Cp (0, 0, 0. 0), Cs = Cs'(0,0. 0, 0). En rappelant que 
QAU J, K Q est Telement [Mi -1) + J. N{k - 1) + Q de la matrice Qx , on en 
deduit: 

15 Q.^AqCsAq'' (23) 

En posant que Qxli j, K Q est Felement lN{i -1) + I , N{k - 1) + jl de Q^on 
obtient la matrice Q^. qui s'ecrit de la fagon suivante : 

= BqCsBq'' (24) 

Btape de blanchiment 

20 La premiere etape du procede cyclique orthonormalise les 

colonnes des matrices Aq ou Bq contenues dans les matrices Qx ou des 
expressions (23) (24). Les matrices Qx et sont hermitiennes de rang P et 
veriflent la decomposition en elements propres suivantes : 

a = et = A, (25) 

25 Ou A ^ est la matrice diagonale de dimension (P x P) des P valeurs non 
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nuUes de et E, est la matrice de dimension [N^xPj des vecteurs propres 
associes. Pour une matrice Bg de rang plein, on peut montrer qu'il y a 
equivalence entecgsupposer que les kurtosis des sources sont de meme 
signe 8 {8 = ±1) et supposer que les valeurs propres de A, sont aussi de 
5 meme signe e. Dans ce contexte. il est possible de construire la matrice de 
blanchlment suivante f de dimension (PxN^): 

T =(A,r£/ (26) 

ou ( A_,) est I'inverse de la racine carre de A^. Des expressions (24) et 
(25) on en deduit que : 
^0 ^TQ, f" = f BQ(sa)Bg"f'' = I, (27) 

Cette demiere expression montre que la matrice f Bg (.eCsV^^ de 

dimension (P x Pi est une matrice unitaire f7 . On en deduit alors que : 

T Bq= U (eQ) (28) 

On rappelle que la matrice Tde blanchiment de Q>, verifle : 

TAq= UXeCS"^ - (29) 

Etape d'isolation cyclique 

On deduit des expressions {28)(29) et (21) que : 

WXo.T,,T2.T3)=TQ,\a,r,.T2.T3)f = U.{eCsfCs\a,Xi,X2,'i{){€Csf if 

WA(X,T,. T2, Ts)=T Q.\a,r,. T2, rs)f^ = U ieCs)"^Cs\a,XuT2,x,) (^C.)"'" £/ ^ ' (30) 
En presence de Pi sources verifiant cf[a,xx,X2,X3)j!tO (1< i < pj la 
matrice (a, rj.r2.T3) est de rang Pi<P. Dans ces conditions les matrices 
xmitaires U et C7 de dimension PxP peuvent se decomposer en deux sous 
matrices de dimension P>cP, et P<(P - Pi) tel que : 

U=[Ux U2] et U =[U, U^] (31) 
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It 

ou les matrices Ui et sont de dimension PkPi et U2 et U2 sont de 
dimension Pk{P-Pi), Les matrices Ui et contiennent les vecteurs 
singuliers associes aux valeurs singulieres non nulles de Wi/ (a, 71,^2,73). On 
en deduit alors que : 

(32) 

Ou la matrice C^^{a.Ti,T2,t3) est une matrice diagonale de dimensions PixPi 
10 compose des elements diagonaux Ci^{a,xi,T2,T3) non nul de Cs^{a,Ti,T2,T3). La 
matrice = (0,0,0,0) de dimension PixPi est composee des elements Ci 
(1< I < Pi) correspondants. Dans ces conditions apres une decomposition en 
valeurs singulieres de Wx^{a,TuT2,T3), on peut determiner les matrices Ti et 7^ 
a partir des vecteurs singuliers associe aux valeurs singulieres non nulles et 
15 T2 et 72 a partir des vecteurs singuliers associe aux valeurs singulieres 
nulles tel que : 

ri=C/i/7i^ , r2=C/2/72", f, = U, /7j"et f,= U,n,'' (33) 

Ou les matrices sont /7i. fh, 11 ^ et TJj sont unitaires. On peut 
construire a partir de Wj^ ia\Ti\T2\T3) une matrice Wf {a\Ti\T2\T3) dependant 
20 uniquement des sources de parametres cyclique (a,xi,X2.T3,e) tel que 
Q''(a,Ti,T2.T3)^0. Pour ce faire on effectue le calcul sulvant : 
f^j (a', ry >,>/)= r," wM\rj\r2\T3) J, = /7, (e Cj"''' C, \aW.T2\Ts) (^Cj"'''y7,» 

J^/(a',r;;r,;T/)= J," W.\a,Tnr2\Ts) f,*= 77, (e C,)"''' C,\a\Tj\r2]Ts) (eC,)''^ n,'' 

JF;(a>/,r,>/)= 7;"jr,^(a',r/>,;r;) f,'=n, {eCS''''c,\a:,Tj\T2,rs){eCj''^n,'' 
25 . (34) 
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De meme on peut construire a paxtir de Wx^'(a\Ti\T2\T3) une 
matrice Wf (a\Ti\r2\T3) ne dependant pas des sources de parametres 
cyclique (a,Ti.T2,T^;e) tel que Ci^(a,Ti,T2.T3)=0 : Autres sources. Pour ce faire on 
effectue le calcul suivant : 

: T2\ r/)= r2" wM\tj : T,: r/) f/ = 7/2 y''^ ^aVr; r/, r, ) 77^ ^ 

(35) 

Ou la matrice C/'(a\ri',r2',r3') est une matrice diagonale de dimensions {P- 
10 Pi)x(P-Pi) compose des elements diagonaux c^ia\Ti\r2\T3') tel que les 
cf[o,,Ti,T2,T3) correspondants sont nuls. La matrice C, de dimension (P- 
Pi)x(P-Pi) est compose des elements o {1< i < Pi) correspondants. 

En particulier dans la T^re version de cette variante on peut 
efifectue Tisolation cyclique en a'=0 et =1. En posant W[ti\t2\t3) = 

15 W/\a\Ti\T2.T3) , C^.(Tj\T2\T3)=C/{a\Tj\T2\T3) Ct CjTi\r2\T3)= C/(a\Tj\T2\T3), 

les equations (34) et (35) devlennent : 

^^A^i\t2\t3)==Ti''wAtj\t2\ts)Ti - 77i Cj"'''c,(ry;r2:irj)(£Cj"''^77," 
(36) 

. ^?;(r/:r/,r/)=r2" fF,(r/:r,;r5)r2 = CJ^''' C,(r/,r/,r3) (f^j"''' 7/2" 

20 (37) 

Etape dldentification 

Les equations (34) et (36) moiitrent que Ton peut identifier les 
matrices unitaires 77i et 77, associees aux sources de parametres cyclique 
25 (a,TuT2,T3,e) en effectuant la SVD{Singular Value Decomposition) conjointe 
des matrices Wf ((xI^ti J,T2J.T3J) pour 1</<K Ainsi pour estimer la matrice 
unitaire de gauche on effectue la diagonalisation conjointe des matrices : 
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Wf (a!M r2W) Wf (cc^.r/^, r/, zs'f pour \<i<K ^38) 

et pour estimer la matrices unitaires de droites on effectue la 
diagonalisatlon conjointe des matrices : - 

Wf {d,rj\T2^,T3Y io!\TjiT2'\T3') pour l<f<K (39) 

5 Pour estimer les matrices unitaires IT2 et associees aux 

sources non associees aux parametres cyclique (a,Ti,T2tT3,e). on eflfectue la 
SVD conjointe des matrices f^^/' (a/, ij J, ^, ^) pour 1<;<K en diagonalisant 
conjolntement les matrices Wf {af^rj J.T2 J.ts J) {aJ,Ti It2J,T3J)^ puis les 
matrices Wf {a^,TiJ,T2J,T3J)^ W^' (a),TiJ,T2^,r3^. 
10 Cormaissant /7i et Ih. on deduit de I'equation (33), a une matrice 

de permutation pres, les matrices unitaires C/i, U2 et U en effectuant : 

Uy^TxHi , f/2= 7^/72 et t/i (40) 

On peut deduire des equations (9) et (29) la matrice Aq a une 
matrice diagonale et de permutation pres. tel que : 

15 f.U^[b,.,,bp\^E,A'''' ^Aq{sCs)'^' AU (41) 

ou est la pseudo-inverse de la matrice T. Chaque colonne , hi (1< I < f5» de 
la matrice 'f U est associe a un vecteur iiq |Cq|l/2 (a^Oq*), 1< q < P. ou fXq est 
un scalaire complexe tel que ||iiq| = 1. En consequence, en transformant 
chaque colonne bi de la matrice 7^ 17 en une matrice Bi de dimension [N x N) 
20 tel que BilU J] = bi[{i -l)Ar + J) (1 ^ t J < iV) on en deduit que: 

Bi-=iig\cg\^^agag^ pour {\<l,q<F) (42) 

Dans ce contexte le vecteur directeur Oq de la cf^^^ source est associe au 
vecteur propre de jB^ associe a la plus forte valeur propre. 

25 Recapitulation de la premiere version du procede cyclique 

En resume, les etapes de la I'^re version du procede cyclique sont 
resumees ci-dessous et sont appliquees sur L observations x(lTe) (I< Z < L) 
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des signaux regus sur les capteurs ( Te: periode d'echantillonnage). 
Estimation 

Etape-1: Estimation des matrices Qx et Q^a partlr des L observations 
x(lTe)* L'estimation de ces, matrices va dependre des hypotheses suivantes : 
5 • Cos stationnaire et centres : Estimateur empirique utilise dans la reference 
[31. 

• Cos cyclo-stationnaire et centres : Estimateur mis en oeuvre dans la 
reference [10]. 

• Cos cyclo-stationnaire et non-centres : Estimateur mis en oeuvre dans la 
10 reference 111]. 

Balnchiment 

Etape'-2 : Decomposition en elements propres des estimees des matrices Qx 
et Qj,. A partir de ces decompositions, estimation du nombre de sources P 
et utilisation des P valeurs propres principales tel que : 
■J 5 O^^ExAx Ex^ - A]^ E^^ ou Ax et A^sont des matrices diagonales 

contenant les P valeurs propres de plus fort module et Ex et E^ sont Igs 
matrices contenant les vecteurs propres associes. 

Etape-3 : Construction des matrices de blanchiment : T = (AJ^'^' E^^^ et 
f = (AJ"'/^ 

20 Etape-4 : Selection des parametres cycliques (a, ti ,T2 ,X3 ,£ ) et estimation de 
la matrice Qx[a,Tut2,T3) a partir des L observations x{lTe). ^estimation de 
cette matrice va dependre des hypotheses suivantes sur les signaux- 
• Cos stationnaire et centres : Estimateur empirique utilise dans la reference 
[3]. 

25 • Cos cyclo-stationnaire et centres : Estimateur mis en oeuvre dans la 
reference [10]. 
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• Cos cyclo-stationnaire et non-centres : Estimateur mis en oeuvre dans la 
reference (11). 

Etape-5 : Calcul d'line matrice Wx^{a,TuT2.T3) de (30ra partir des matrices 
QxHoi,Tur2,T3) . Tet r. Apres une decomposition en valeurs singulieres de 
5 Wj{a,Ti,T2sT3] , determination des matrices unitaires Ti et T, associe aux 
valeurs singulieres non nulles et T2 et Tj associe aux valeurs singulieres 
nuUes. 
Selection 

£tape-6 : Selection de K triplets de retards (xi^T2^T3'^) ou 

10 |Ti'f| + |T2'^| + |T3'^|9tO. 

Estimation 

Etape-7: Estimation des K matrices Qa<t^i^T2^,X3^) de (2). Comme dans 
Tetape-l cette estimation va dependre des hypotheses faite sur le signal tel 
que : 

15 • Cos stationnatre et centres : Estimateur empirique utilise dans [3]. 

• Cos cyclo-stationnaire et centres : Estimateur mis en oeuvre dans [10]. 

• Cos cyclo-stationnatre et non-centres : Estimateur mis en oeuvre dans [1 1], 
Identification 

Etape-8 : Calcul des matrices Ti 0x(Ti*^»T2^T3'^) Ti» et estimation de la 
20 matrice unitaire Ui (associe aux parametres cycliques (a, xi ,X2 ,T3 en 
effectuant la diagonalisation conjointe des K matrices : 

Etape-9 : Calcul des matrices T2 Qxi'^i^. '12^,^3^] T2" et estimation de la 
matrice unitaire U2 (associe aux autres sources) en effectuant la 
25 diagonalisation conjointe des K matrices T2 Qx('i^i\T2'',T3'') T2". 

Etape-lO : Calcul de la matrice imitaire U en effectuant : i7=[L7i U2] 
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Etape-11: Calcul de T^U =[bi...bp] et construction des matrices Bi de 
dimension {2V x JV) a partir des colonnes bi de T*l/. 

Btape-12: Estim^gpn des signatures a<j {!< q < F) des P sources en 
appliquant une decomposition en element sur chaque matrice Bi 

Recapitulation de la deuxiezne version du procede cyclique 

Les etapes de la 2^^^^ version de la methode FOBIUM cyclique sont 
restimees ci-dessous et sont appliquees sur L observations x{lTe) (1< IS L) 
des signatix regus sur les capteurs ( T^: periode d'echantillonnage). 

Etape-l: Estimation des estimees des matrices Qx et O^a. partir des L 
observations x(lTe). Uestimation de ces matrices va dependre des 

- a*- 

hypotheses suivantes : 

• Cos stationnatre et centres : Estimateur empirique utilise dans [3]. 

• Cos cyclo-stationnaire et centres : Estimateur mis en oeuvre dans [10], 

• Cos cyclo-stationnaire et non-centres : Estimateur mis en oeuvre dans [1 IJ. 
£taper2: Decomposition en elements propres des matrices Qx et Q^. A 
partir de ces decompositions, estimation du nombre de sources P et 
utilisation des P valeurs propres principales tel que : Q^Ex Ax E^^ et = 

" ou Ax et sont des matrices diagonales contenant les P valeurs 

propres de plus fort module et Ex et £^sont les^' matrices contenant les 

vecteurs propres associes. , V 

Etape-3: Constmction des matrices de blanchiment : T = (AjJ^'^Ec" et 
f=(A,r/'£^H 

Etape'4: Selection des parametres cycliques (a, xi ,T2 .xs ,e ) et estimation de 
la matrice Qx[a.Tj,T2,T3) a partir des L observations x{iTe), L'estimation de 
cette matrice va dependre des hypotheses suivantes sur les signaux : 
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• Cos stationnatre et centres : Estimateur empirique utilise dans [SJ. 

• Cos cyclo-stationnaire et centres : Estimateur mis en oeuvre dans [10], 

• Cos cyclo-stationnaire et non-centres : Estimateur-mis en oeuvre dans [1 1]. 
Etape-5: Calcul d'une matrice Wx'(a.rj, 72,^3) de (30) a partir des matrices 

5 Qx[a,TuT2.T3] , Tet r. Apres une decomposition en valeurs singulieres de 
W/(a,ri,r2.r3) , determination des matrices unitalres Ti et associe aux 
valeurs singulieres non nulles et T2 et 7^ associe aux valeurs singulieres 
nulles. 

£tape-6: Selection de K ensembles de parametres {a^Ti^T2^T3^e'^). 
10 £tape-7: Estimation des K matrices Qjt\aKxiKT2Kx^^) de (19). Comme dans 
Tetape-l cette estimation va dependre des hypotheses faite sur le signal tel 
que : 

• Cos stationnatre et centres : Estimateur empirique utilise dans [3]. 

• Cos cyclo-stationnaire et centres : Estimateur mis en oeuvre dans [lOJ. 

15 • Cos cyclo-stationnaire et non-centres : Estimateur mis en oeuvre dans [11], 
Etape-8 : Calcul des matrices Ti g/^(a ^xi^,T2^T3'^) TiH et estimation de la 
matrice unitaire Ui ou f/, {associe aux parametres cycliques (a. xi ,T2 .T3 ,8 )) 
en effectuant la SVD conjointe des K matrices : Ti Qx^H(xK^iKr2Kx3^) Ti", 
Etape-9 : Calcul des matrices T2 Qx^Ha ^xi^, T2^T3^-) et estimation de la 

20 matrice unitaire U2 ou (associe aux parametres cycliques (a, xi .X2 ,X3 ,8 )) 
en effectuant la SVD conjointe des K matrices : T2 P/*'(a^Xi^X2^T3fc) T2". 
Etape-lO : Calcul de la matrice imitaire C/en effectuant : U=lUi U2] 
Etape-11: Calcul de T^U =fbi...bp] et construction des matrices Bi de 
dimension (iV x iV) a partir des colonnes bi de'fu. 
25 Etape-12: Estimation des signatures Oq (1< q < ^es P sources en 
appliquant une decomposition en element sur chaque matrice Bi. 
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REVENDI6ATIONS 

1 - Procede dldentlficatlon autodidacte a Tordre 4 de signatures d'au moins 
deux sources dans un systeme comportant un nombre de sources P et un 
nombre N de capteurs de reception recevant les observations, lesdites 
sources ayant des tri-spectres differents caracterise en ce qu*il comporte au 
moins les etapes suivantes : 

a) une etape de blanchiment a Tordre 4 des observations re?ues sur les 
capteurs de reception afin d*orthonormaliser les vecteurs directeurs des 
sources dans les matrices de quadricovariance des observations 
exploitees, 

b) une etape de dlagonalisation conjointe de plusleurs matrices de 
quadricovariance blanchies (etape a) afin d'identifier les signatures 
spatiales des sources. 

2 - Procede selon la revendication 1 caracterise en ce que les observations 
utilisees correspondent aux matrices de quadricovariance moyennees 
temporellement definies par 

p 

Qx{iui2,X:;) = ^y^(^i.'^2,T3) {ap® Op *) {Op ® Up y (4a) 

= Aq G(Ti,T2,T3) Aq^ (4b) 

ou Aq est une matrice de dimension (iV^x P) definie par Aq = [{ai®ai*), 

(ap®ap')I, Cs(ti,T2,T3) est une matrice diagonale de dimension (P x F9 deflxiie 
par Cs(xi,T2,T3) = diagIci(Ti,X2»T3)»..., Cp(ti,T2.T3)] et ou Cp(ti,T2,T3) est defini par : 

Cp(Ti ,T2,T3) = <Cum(jr^(0, Sp{t-Xx)\ Spit-Xi), Sp(t-T3))> (5) 

3 - Procede selon la revendication 2 caracterise en ce qu'il comporte au 
moins les etapes suivantes 

Etape 1 : estimer, par ^x, la matrice Qx. a partir des L observations x{lTe) 
en utilisant un estimateur non biaise et as3rmptotiquement consistant. 
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Btape 2 : decomposer en elements propres ^x, estimer le nombre de 
sources P et restreindre la decomposition en elements propres aux P 
composantes principales : « fex fL t^, ou Ax est la matrice diagonale 
contenant les P valeurs propres de plus fort module et tx est la matrice 
5 contenant les vecteurs propres associes. 

Etape 3 : construlre la matrice de blanchlment : t = (A^)'^^' 

Btape 4 : selectionner K triplets de retards {xiKx-zKx^f^) ou 

|Tl'^[ + |T2'^| + |T3fci9£0. 

Etape 5 : estimer , par ^x(xi^T2^,T3'^). les K matrices Qo<k'^iKi2Kxz% 
10 Etape 6 : calculer les matrices t ^JixiKx^Kxz^) ^ et estimer, tJs6L. la matrice 
unltaire Usoi par diagonalisation conjointe des K matrices ^ ^x(xi^T2^T3'c) t« 
Etape 7 : calculer t *£/soz=[6i,,.ftp] et construire des matrices hi de 
dimension [NkN). 

Etape 8: estimer. par hp, les signatures Og (1< q < P) des P sources en 
15 appliquant une decomposition en element sur chaque matrice 6l 

4 - Procede selon I'une des revendications 1 a 3 caracterise en ce qu*il 
comporte au moins une etape d'evaluation de la quaHte de Tidentiflcation 
du vecteur directeur associe en utillsant un critere tel que 
20 D{A A) = (ai. a2 , Up) (16) 

ou 

ttp = mm[d(ap, d<)] nyi 

l<i<P V ■»> « i 

et ou d(u,p) est la pseudo-distance entre les vecteurs u et v. tel que : 

^^'^ «») = 1 - 73CiriO (18) 
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5 - Procede selon Tune des revendications 1 a 3 caracteiise en ce qu*il 
comporte au moins une etape desolation cyclique a Tordre 4 apres I'etape a) 
de blanchiment a Tordre 4. 

5 6 - Procede selon la revendication 5 caracterise en ce que Tetape 
d identification est realisee en utilisant des statistiques cycliques d'ordre 4. 

7 - Procede selon Tune des revendications 1 a 6 caracterise en ce que le 
nombre de sources P est superieur ou egal au nombre de capteiu-s. 

10 

8 - Procede selon Tune des revendications 1 a 7 caracterise en ce qu*il 
comporte au moins une etape de goniometrie a parUr de la signature 
identiflee des sources. 

15 9 - Procede selon Tune des revendications 1 a 7 caracterise en ce qu'il 
comporte au moins une etape de flltrage spatial apres la signature identifiee 
des sources- 

10 - Utilisation du procede selon Tune des revendications 1 a 9 a un reseau 
20 de communication. 
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